UNIFORMIZING DIFFERENTIAL EQUATIONS OF ARITHMETIC 

(l;e)-GROUPS 

JORG HOFMANN 

Abstract. We present a numerical method to compute accessory parameters of the uni- 
formizing differential equations of elliptic curves associated to arithmetic Fuchsian (l;e) 
groups. Up to PSL2{M.) conjugation there are only finitely many such groups and we give a 
complete list of the accessory parameters. 



Introduction 

A discrete subgroup T of P5L2(M) is called Fuchsian group and its elements act on the 
upper half plane EI as Mobius transformations. The quotient Y{T) = M/T can always be 
equipped with a complex structure and is therefore a Riemann surface. An important class of 
Fuchsian groups are congruence subgroups of 5L2(Z) which lead to modular curves, which can 
be compactified by adding finitely many cusps. An immediate generalization of congruence 
subgroups are arithmetic Fuchsian groups. A group is called arithmetic if it is commensurable 
with the embedding of the norm one elements of an order of a quaternion algebra into M2 (M) . 
If r is a an arithmetic Fuchsian group, then X{T) can be realized as a projective algebraic 
curve defined by equations with coefficients in a number field. In section [l] we briefly review 
the basic facts on arithmetic Fuchsian groups. Special attention will be paid on the case 
where F is a Fuchsian group of (1; e)-type, that is a Fuchsian group whose quotient has genus 
one and the projection : H — )■ H/F branches exactly at one point and the branching index 
at this point equals e. Arithmetic Fuchsian groups of (l;e) type were classified by Takeuchi 
in |Tak83] . The main issue will be the inverse w of cj). As will be explained in section [2| lo 
can be realized as the quotient of two linearly independent solutions of a differential equation 
L, called the uniformizing differential equation. Recently Sijsling gave a list of equations 
for projective curves that yield models for most of the Riemann surfaces associated to the 
groups from Takeuchis list [SijbJ . If such a model is known, it determines the local data 
of the uniformizing equation. But in general, one complex parameter, called the accessory 
parameter, remains undetermined. In section [3] we explain how a method that computes 
high precision approximations of the generators of the monodromy group of a linear Fuchsian 
differential equation can be used to approach the determination of the accessory parameter. 
In all cases we are able to compute floating point approximations of the accessory parameter 
with accuracy sufficiently high to identify the accessory parameters as algebraic numbers. 



Date: March 4, 2013. 

1991 Mathematics Subject Classification. 34M35, 33E10 ,20H10. 

Key words and phrases. Lame equations, arithmetic Fuchsian groups, Uniformization. 



UNIFORMIZING DIFFERENTIAL EQUATIONS OF ARITHMETIC (l;e)-GROUPS 



2 



1. RiEMANN Surfaces and Arithmetic Fuchsian groups of (l;e) type 

One way to construct Riemann surfaces is by taking the quotient M/T by a suitable group 
acting on the upper half plane. We will focus on the case where T is an arithmetic Fuchsian 
group of (1; e)-type. The books |Vig80| and |Kat92j cover all the material needed. 
We start by recalling the necessary definitions. For a Fuchsian group F whose action on IHI 
has a compact fundamental domain the quotient M/T is also compact with genus g. Since 
then necessarily F has no parabolic elements, if the projection H/F — )• IHI branches above r 
points with index mi, . . . , rrir the signature of F can be defined as the tupel {g; nii, . . . , nir). 
Groups with signature (0; mi, m2, 711.3) are called triangle groups and their investigation has 
a long history. Beyond triangle groups a prominent role is played by groups with signature 
(l;e). 

Definition 1.1 (Fuchsian (1; e)-group). A subgroup F < PSL2{^) is called Fuchsian (l;e)- 
group if the quotient H/F of the action of F on the complex upper half plane is a compact 
Riemann surface of genus one and the projection HI — )• H/F ramifies of index e above exactly 
one point. 

Following Fricke and Klein |FK65j . a lift of a Fuchsian (l;e)-group F to S'L2(M) can be 
presented as 

< a, /3, 7|a/3Q!~"^/3~"^7 = —id2,j'^ = —id2 >, 

where a and /? are hyperbolic elements of S'L2(M) and 7 is elliptic with tr(-y) = 2cos(^). In 
addition the Fricke relation 

tr{af + tr{pf + tr{al3f - tr{a)tr{f3)tr{a/3) = 2 - 2cos (J^^ 

holds. For realizations A,B G SL2M. of a and /3 and p the fixed point of the elliptic element 
ABA~^B~^, the matrices A and B considered as maps are the side pairings of the hyperbolic 
polygon P with vertices p, Ap, Bp, ABp and angle sum 2^. It is Poincare's theorem which in 
this case states that P is a fundamental domain for F. After suitable choices the fundamental 
domain of a Fuchsian (1; e) group in the upper half plane can be depicted as a quadrilateral, 
as in figure [TJ There are even up to conjugation infinitely many such groups, but if we impose 




Figure 1. Fundamental domain of a (l;e) group 

the condition that the group is arithmetic, a notion we will make precise below, only finitely 
many conjugation classes exist. The notion of quaternion algebras is crucial. 

Definition 1.2 (quaternion-algebra). For two elements a,b of a field F of characteristic not 
equal to 2, the algebra Fi © Fj © Fk © Fl with multiplication given by i'^ = a, = b and 
ij = ji = k is denoted by (^) and called the quaternion algebra determined by a and b. 
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Two familiar examples of quaternion algebras are the Hamiltonian quaternion algebra 
and the matrix algebra (^) ~ M2{F). The next lemma shows that after an extension 
to the real numbers, quaternion algebras are governed by the above two examples. 

Lemma 1.3. Let (pi be an embedding of the number field F in M, then ^F,4ii is either 
isomorphic to M2(M) or a division algebra. 

If j M = M2(M) the embedding 0j is named split and nonsplit elsewise. In the 

following let F be a totally real number field of degree n. There are n distinct embeddings 
(pi of F into the real numbers, called infinite places. Consider quaternion algebras which are 
split at the identity (pi and non split at all other infinite places. Denote the corresponding 
isomorphism M — ?■ M2(M) by j. 

Definition 1.4 (Order). The ring of integers of F are denoted Zp. We say that a Zp 
submodule of 'S^F,i ^ is an order if it is also a subring. 

Clearly an element x of a quaternion algebra (^) can be written as x = foi + fij + f2k + fsl 
and 

\x\ :=f^-aff-bfl+abfi 
introduces a norm. Denote all elements of norm 1 of an order O by O^. 

Definition 1.5 (arithmetic Fuchsian group). A Fuchsian (l;e)-group F is called arithmetic 
if there exists a quaternion algebra with order O such that F is commensurable with j(O^). 

Takeuchi used that an arithmetic Fuchsian (1; e)-group F can be described completely in 
terms of the triple 

{tr{a),tr{f3),tr{al3)) 

to prove that up to i-*S'L2(M)-conjugation there are 73 arithmetic Fuchsian groups of (l;e)- 
type. A set of generators (oq, Pq, 70) of F is called fundamental, if 

2 < tr{ao) < tr{l3o) < tr(ao/3o) 
and if it is of least height among all such generating triples, where the height is defined as 

h{a, /3, 7) := tr{af + tr{l3^) + tr{a(3f. 

In |Tak83| a list of these fundamental triples {x,y,z) for arithmetic Fuchsian (l;e)-groups is 
given. Moreover, the quotient H/F is an Riemann surface with a marked point, the point 
above which the projection EI — )• H/F branches. Hence this quotient is an instance of an 
orbifold. 

2. Orbifold Uniformization and Differential Equations 

In this section we recall shortly the classical theory of orbifold uniformization. A reference 
which explains the mathematical as well as the historical aspect of the uniformization of 
Riemann surfaces is the book |dSG| by the collective Henri Paul Saint Gervais. 

Definition 2.1 (orbifold). Let X be a complex manifold and Y C X a hypersurface, that 
splits as y = DjYj in irreducible components. Furthermore associate to every Yj a natural 
number bj > 2 or 00. The triple {X,Y, (bj)j) is called an orbifold if for every point in 
X\L)j {Yj\bj = 00} there is an open neighbourhood U and a covering manifold which ramifies 
along U nY with branching indices given by {bj)j. 
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If the above local coverings can be realized globally this definition can be subsumed as. 

Definition 2.2 (uniformization). If there exist a complex manifold M and a map : M — )• X 
which ramifies exactly along the hypersurfaces Yj with the given braching indices hj , then the 
orbifold {X,Y, is called uniformizable and M is called a uniformization. 

The multivalued inverse of (j) is called developing map and it is in general a hard problem 
to determine this map explicitly. We will focus on the case where X is the quotient of the 
upper half plane by an arithmetic Fuchsian group of (l;e)-type and therefore M coincides 
with upper half plane itself. The hypersurface Y is just one point above which the map <j) 
ramifies with index e. There is a fundamental connection between the multivalued inverse of 
the map <j) and a certain linear differential equation, which goes back to Poincare. For a local 
coordinate x on X the inverse map u = (j)~^{x) is PGL2(M)-multivalued, but the Schwarzian 
derivative 

/ /// \ o / " \ ^ 



{Sw) :-- 



is a single- valued map on X. For a proof of the last two claims and the subsequent desired 
proposition connecting the Schwarzian derivative and differential equations see |Yos87j [Ch.4.2] . 

Proposition 2.3. Let oj{x) he a non-constant PGL2(^) -multivalued map, then there exist 
two C-linearly independent solutions yo{x),yi{x) of the differential equation 

d^ 1 1 

+ -.S{io{x)) 



such that Loix) = . 

Definition 2.4 (uniformizing differential equation). The differential equation from the above 
proposition is called uniformizing differential equation or Schwarzian differential equation. 

In general the unformizing differential equation can not be computed from F directly, but 
one can try to recover Le from its monodromy group M, since the projectivization of M 
coincides with F. To be more precise we recall the definition of the monodromy group of a 
linear homogeneous differential equation L. Such an equation L of order n on a Riemann 
surface S has near any ordinary or regular singular point p exactly n over C linear independent 
solutions, denote the space spanned by these solutions by Fp. Since analytic continuation of a 
germ of Os along a path on S is defined up to homotopy, there is a well defined representation 

p:7ri(S\S,p)^GL„(C), 

where S is the set of singular points of L. For a loop 7 G 5* with base point p, analytic 
continuation along 7 transforms Fp into Fp and since Fp and Fp are two bases of the same 
vector space, the solution space of L, analytic continuation yields an element M G G/„(C) as 

Fp = FpM. 

Definition 2.5 (monodromy representation, monodromy group). The representation p is 
called monodromy representation of L and its image in GLn{C) is called the monodromy 
group. 

Denote the set of branch points of the projection M ^ X hy B and let p a point outside 
of B, the representation pQ : 7ri{X \ B,p) — )• F C PSL2{^) associated to the covering 
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H — 7- X is conjugated to the projectivization of the monodromy representation of p : 'Ki{X \ 
B,p) — )• GL2(C) — PGL2{C). This can be seen as follows. Two solutions yi,yo of Le are 
a local branch of the inverse u}{x) of the covering map. By analytic continuation along a 
loop 7 centred at p two C-linear independent local solutions yi , yo oi Le are converted to 
m = ayi + byo, yi = cyi + dyo with M = ^) G GL2(C). This yields 

yi ^ ayi/yQ + h 
yo cyi/yo + d' 

which is again the quotient of two solutions of Le and therefore another branch of ij{x) related 
to the former one by the Mobius transformation associated to M. Hence if the monodromy 
representation is computed according to an arbitrary basis Fq the two representations po and 
p are conjugated. 

We will do no computations with L e directly, instead we consider a differential equation Lpi 
on the projective line which is related to via a certain pull-back. Since we assumed X to 
be a compact Riemann surface it can be realized as an elliptic curve, that in an affine chart 
can be given in the form 

E : y'^ = 4x^ + ax + b, a,b e C. 

The map tt :£'—;■ \ S, (x, y) i— )• x is a twofold cover of \ {xi, X2, X3, 00}, branching of 
index two at the points of S := {xi, X2, X3, 00}, where Xi are the roots of P(x) := 4x^ + ax + 6. 
Thus Le IS the pullback of a differential equation on P^, furnished with a local coordinate t 
which is defined by local data as 

1 1 



Lm 



y = 0, n 

^ ' 2e 2 



This differential equations is called algebraic Lame equation and the complex number C is 
referred to as accessory parameter. This kind of differential equation was first investigated 
in connection with ellipsoidal harmonics |Poo07j [ch. IX] . The accessory parameter does not 
affect the local exponents and the location of the singularities of Lpi . That Lpi is indeed the 
correct object can be explained by a closer look at its Riemann scheme 
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A basis of the solution space of Lpi at a regular point can be given by two power series and 
in a vicinity of any of the finite singular points it is given by 



yo(i) = - ij)" and 

n=0 

00 

yi{t) = it-ti)^^bn{t-ti 



n=0 



whose quotient is = {t -ti)^ X]^o'^"(^ ~ ^0)" CHcMj. Similarly since -§ - = ~2i 
the shape of this quotient at infinity is 
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At an ordinary point p of Lpi the local exponents are and 1, hence following |Inc56j [p.396 
ff] in a vicinity of p the Ansatz 

oo 

y^{t) = ^ an{a){t - pr+'' = 0, a e {0, 1} 

provides a basis of the solution space of Lpi consisting of convergent power series, adding 
derivatives a fundamental matrix can be found as 



The radius of convergence equals the minimal distance from p to any of the singular 

points of Lpi . The coefficients are computed recursively as 

52((i - 1 + ct)0- - 2 + a) + (j - 1 + a))aj.i 



a„((T) 



g3{n + a){n - 1 + a) 
(n(n + 1) + C - 4(i - 2 + - 3 + a))aj-2 
gs^n + a){n -l + a) 
(6(j - 3 + a))a,- 



n > 3. 



fl'3(ra + cr){n -l + a) 

where the initial conditions, i.e. the values of oq, ai and 02 could be chosen arbitrarily with at 
least one Oj 7^ 0, but we will fix them as oq = 1, ai = and 02 = 0. The monodromy group of 
Lpi is generated by the images of the monodromy presentation of loops 71,72,73 starting at 
p and encircling exactly one of the finite points in S in counterclockwise direction, with the 
additional property that composition of pathes 717273 is homotopic to a path 700 encircling 
00 once in clockwise direction. The images of theses pathes under p are called Mi,M2,M3 
and Moo- To compute the analytic continuation choose points po, • • • ,Pm-i,Pm = Po on (a 
path homotopic to) 7 s.t. pi+i lies in the radius of convergence of the elements of Fp.. With 
the above choice of initial conditions Fp-(pi) = id2 holds and 

is implied. Which finally leads to 

m—l 

= Yl Fp„_,(Pm-i-l). 
i=0 

To relate the monodromy representation of Lpi with our initial group T it is enough to 
understand how the monodromy group of \S lifts to the fundamental group of X. We recall 
the standard construction of the fundamental group of an pointed elliptic curve as branched 
two fold covering of the four times punctured sphere. Cut the Riemann sphere from ti to t2 and 
from ^2 to c«, expand the cuts a little and glue two copies of this cutted sphere along the cuts 
with opposite orientation. The fundamental group of the pointed elliptic curve is generated 
by lifts of the loops 5 and 7 as in figure [2| Hence the monodromy of Le along 5 coincides 
with the monodromy of Lpi along a loop 7172 and the monodromy oi Le along 7 coincides 
with the monodromy of Lpi along 7oo72^^ ~ 73 ° 7i- Thus we can identify up to scalars F 
with the group generated by M3M1 and M1M2. Note that M3M1M1M2 = M3M2. According 
to Fricke the traces of Mj, i = 1, . . . , 3 and the traces of the products above determine the 
monodromy representation of Lpi. This relation between the monodromy of Lpi and the 
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Figure 2. Connection of tti{E) and 7ri(P^ \ S) 



Fuchsian group T can also been established, in the spirit of Klein and Poincare, by a look at 
the image of a quotient y = of the upper half. Then suitable analytic continuation of y and 
the Schwarz reflection principle would exhibit how certain elements of the monodromy group 
of Lpi coincide with side pairings of the fundamental domain of F. This is exactly the way 
how the monodromy of the hypergeometric function is related to triangle groups. Four points 
on the Riemann sphere can always be mapped to 0, oo by a Mobius transformation, 
that means given a Fuchsian (1; e)-group Lpi depends on two parameters, namely A and the 
accessory parameter. It is explained in |Iha74| that A and C have to be elements of Q and 
moreover in a suitable coordinate they have to be elements of a number field that can be 
determined from F directly. But recently J. Sijsling |Sijb| constructed WeierstraB-equations 
for members of isomorphism or at least isogeny classes of the elliptic curves associated to 
those 73 groups which provides the local data of the uniformizing differential equation i.e. the 
location of the singular points and the Riemann scheme of the corresponding Lame equations 
are known in these cases. Whenever F is commensurable with a triangle group he additionally 
used in [Sija] Belyi maps to determine the accessory parameters of the uniformizing differential 
equations. We will use his results and numerical methods in the next section to tackle the 
problem of the determination of accessory parameters by investigating the monodromy of the 
uniformizing differential equations. 



3. Approximation of Accessory Parameters 



In |CC88] G.V. and D.V. Chudnovsky considered algorithms to compute approximations 
of generators of the monodromy group of linear differential equations. They also had avail- 
able implementations and used them to do high precision approximations to guess accessory 
parameters in at least a few cases picked from the Lame equations associated to Takeuchis 
list |CC89| . In a talk in Banff in 2010 [BeulOj F. Beukers stated that Sijsling used similar 
numerical methods and was again successful in a few cases. But to the authors knowledge 
there is no complete list of accessory parameters available. Hence we picked up that approach 
and did our own implementations as follows. Once we gained the insights from the last section 
the approach is straight forward. Replace the solutions yi and yo of 



Lml 



y = 0, n 



1 

2e 
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by truncations 

N 

yf(t) = ^a„(a)(t-pr+'^ = 0, a e {0,1} 

i=0 

and the fundamental matrix Fp by 

and approximate the monodromy along the path 7 as 

m— 1 
i=0 

How this approximation is done effectively was investigated by G.V and D.V. Chudnovskys in 
|CC87| and |CC88) . There they considered and partially answered questions dealing with fast 
computations of the coefficients a„(cj) and how the points pi should be chosen. For details 
on the choices of loops representing a homotopy class and vertices of the involved polygons 
and only the implementations done see |Hof| . The errors occurring due to the truncation 
process can be bounded by explicit bounds. As explained above the monodromy group of 
is generated by the products M1M2 and M3M2. In most cases the local data of Lpi is 
known by Sij slings work, and the determination of the accessory parameter is left. It has to 
be determined such that the trace triple 

(1) {tr{MsM2),tr{M,M2),tr{M3Mi)) 

or a permutation hereof coincides with one of the triples in Takeuchis list. If only a permu- 
tation of the trace triple is found, one can use a Mobius transformation interchanging the 
singular points of Lpi to recover the actual triple. If only the isogeny class of E is known, we 
started with one member of this class and computed all elliptic curves E, such that there is an 
isogeny p : E ^ E oi degree bounded by a given n € N and tried to find algebraic accessory 
parameters for differential equations with the local data specified by E. The tools provided 
by |Sijc| proved to be very helpful. Denote the tuple of generators of the monodromy of Lpi 
with accessory parameter C by M{C) = {Mi{C) , M2{C) , M^{C)) and its approximations by 
A{C) = {Ai{C) , A2{C) , A'iiC)) . Since the monodromy group for the correct value of C is a 
subgroup of SL2{^), it is reasonable to consider the maps 

/(T) ■.= max{Im{tr{TiT2)),Im{tr{T?,Ti)),Im{tr{T^T2))] and 
g{T) := \\Im{tr{TiT2))\\ + ||/m(tr(r3Ti))|| + ||/m(tr(r3r2))|| 

for a 3-tuple T of GL2{C) matrices. By a theorem of Hilb |IIil08| . there is an infinite but 
discrete subset of accessory parameters in the complex plane such that /(M(C)) = and 
g{M{C)) = 0, but luckily numerical investigations suggest, that only the one which has 
smallest absolute value is algebraic. If the coefficients of the elliptic curve E are real, we use 
the function / to approximate C. All computations are done in multiprecision arithmetic as 
provided by the C library mpc [EGTZll] or the computer algebra system maple. Start by 
choosing a small number ex and Ci, such that the signs of f{N{0)) and f{N{Ci)) differ and 
build the pair S = (SI, 52) = {f{N{0)),f{N{Ci))) which is used as input of algorithm [ij 
If the accessory parameter is expected not to be real, use the map g. Choose a small real 
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Algorithm 1: Algorithm for real coefficients 
input: S 

if \ f{N{Ci)\ < ex then 
I return Ci 

else if |/(iV(0)| < ex then 
I return 
else 

t = Ci/2 

while \f{N{t)\ > ex do 
t = (51 + S2)/2 

if sgnifiN{t))) = sgn{f{N{Sl))) then 

S=(t,S2) 
else 

S=(Sl,t) 
end 

adjust N and the length of the mantissa 
end 

return t 
end 



number r and C = to compute g{N{C)) and continue with the code from algorithm [2| 
Adjust N and the mantissa means that in every step the number of coefficients in the power 



Algorithm 2: Algorithm for complex coefficients 
input: t,C 

if g{N{C)) < then 

I return C; 
else 

while g{N{C)) > ex do 
cont = false; 
for /c = to 7 do 

CC = C + rexp{2iTi/k) 
if g{N{CC) < 5(A(C)) then 
C = CC; 
cont = true; 

adjust N and the length of the mantissa 
break; 
end 

if cont = false then 

r = r/2; 
end 
end 
end 

return C; 
end 
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series expansion of the solutions and the length of the mantissa are chosen as short as pos- 
sible to guarantee that the errors made by the truncation and by cutting off the involved 
floating point numbers is smaller than the precision we have already achieved in the traces 
of A{C) compared to the target traces M(C). The precision one needs to identify a floating 
point number as an algebraic number depends on the height and the degree of its minimal 
polynomial. An algorithms that accomplishes the identification of floating point numbers 
as algebraic expressions is described in [FBA96j and there exist various implementations for 
example in maple. Therefore the running time of the algorithms [T] and [2] for given ex depends 
mainly on two factors. Namely on the height of the accessory parameter and in addition on 
the location of the singular points of Lpi . The influence of the first factor is obvious. The 
location of the singular points affects the running time as the configuration of the singular 
points mainly determines the number of expansion points pi needed in the approximation of 
the monodromy as explained above. Hence in the simplest cases after only a view minutes 
we obtained a promising candidate for the accessory parameter, but in the more complex 
cases the computer had to run for several hours. Once we obtained such a candidate Caig we 
checked that A{Caig) and the trace triple from Takeuchis list coincid up to permutation up 
to at least 200 Digits. The reason why we can not turn our method in a rigorouse proof is 
that we do not have any a priori knowledge on the heights of the coefficients of Lpi. 

The elliptic curves and accessory parameters found in this way are listed in the tables 
below. The labels are chosen according to |Sijb| . The symbol 

rid/ nor 

encodes 

• rid', the discriminant of the coefficient field F of (^) 

• ud- the norm of the discriminant of the associated (^) 

• r: roman number used to distinguish cases with equal rid and n/j. 

The data in the tables below is the label nd/nof an equation for the corresponding ellitic 
curve E : + x'^ + aixy + oay = x^ + a2x'^ + a/^x + ag, the accessory parameter C of the 
pulled back uniformizing differential equation Lpi and the the squares of the entries of the 
trace triple ([T]) . The polynomial P which specifies the coefficients of the differential equation 
-Lpi can be recovered from E by the substitution of y by \{y — aix — a^) and the elimination 
of the x^ term. If an algebraic number 7 is involved in the coefficients of E or the given 
trace triples has a complicated radical expression, we give some digits of its fioating point 
expansion and its minimal polynomial Whenever all nonzero singular points of Lpi have 
non- vanishing real part the basepoint p is chosen as the imaginary unit i, elsewise it is chosen 
as 0. The singular points ti, i = 1, . . . ,3, which are the roots of P are ordered according to 
the argument of the complex number ti — p. This fixes the loops 7^ and hence the trace triple 

0- 

Table 1: ramification index 2 



2/2 + xy + y = - 334x - 2368 
n — 79 

'-^ - ~64 

(5,12,15) 
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iP' = — T?' — Ax + 4 

^ = h> 

(8,6,12) 


1/14 


y2 = ^3 1 23220a; - 2285712 

^ - 8 

(7,7,9) 


5/4i 


y2 ^ 3 _ f.l^ , 1W2 _ , 3^W_ 25 _ 11/5 

(3 + ^/5,6 + 2^5, 7 + 3^5) 


5/4ii 


y2 , , 1 y = ^3 _ C_iy5 , 1) ^2 _ C_329 , m /gX ^ , 769 _ 287/5 

^ = "lis + 

(^ + 1^5,6 + 2^5,1 + 1^5) 


5/4iii 


y2 , , 3 _ 1^5 + \\ a;2 _ (101 , 45 ^ _ 1895 _ 847/5 

r; - --^ - i^/5 
~ 64 32 ^0 

(i + iV5,i + iV5,i + iV5) 


8/7i 


y2 = ^3 , ^2:^ - f 142 a/2 + 202') X - 1170a/2 - 1655 

C = -ii-iV2 

(3 + v^, 12 + 8\/2, 13 + 9\/2) 


8/7ii 


y2 = 3,3 _|_ ^/2x2 _ f 142a/2 + 202") x - 1170a/2 - 1655 
(3 + 2\/2, 5 + 3V2, 5 + 3\/2) 


8/2 


2/ = X — X 

C = 

(3 + 2 V2, 5 + 3^2, 5 + 3^2) 


12/3 


y2 = ^3 _ _ 1) ^2 _ (65^ _|_ iii^ rj. _|_ 348v^ + 603 

C7=| + ^V3 

(8 + 4^3, 3 + V3, 9 + 5^3) 


12/2 


y2 = 3,3 ^ ^^2 _^ ^ _,_ 3 ^ _ 5 

C = 

(4 + 2^3, 4 + 2^3, 4 + 2v^) 
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13/36 


y2 + (^3+713^ y^^3 5601845 60077-16383713^ ^ + 1551027^ 

^ 1625 1 375 /To 
^ 128 128 ^ 

(1 + i^/l3, 16 + 4^13, f + |^/l3) 


13/4 


y2 + (^3+713^ y ^ ^3 (^3+713 j ^2 (^275+757l3^ ^ 1565-433713 

65 15 ./To 

/ 11+3713 11+3713 7+7l3\ 
\ 2 ^ 2 ' 2 y 


17/2i 


y2 + + (7 + 1) y = + ax^ - (-6I7 + 157) x + 3487 - 896 

' 2 
^ -55+207 

(^5+7l7^10 + 2v^, 21+5v^) 


17/2ii 


2/^ + xy + (a + 1) y = x^ + ax^ - (-6I7 + 157) x + 3487 - 896 

_ 1+717 

' 2 

^ -55+207 

- 64 

(^+f^,5 + yT7, 13+37I7\ 


21/4 


y2 ^ ^l^+19y+15^ y _ ^3 ^7^-97+15^^2 ^9973-14472-6187-1451 j ^ 

887^-50572-7777-555 
14 

= x^ - 4x=^ - x^ + lOx + 43, 7 = -1.292 - 1.323? 

15-5A/2T+56i73-24i77 
^ 128 

/ I5+372T 5+7773 17+37773\ 
\ 2 ' 2 ' 2 / 


24/3 


y2 = a.3 ( ^1,^3 1^2 l^ + l)a;2 (17295^3 14243^2 604597 + 155218) X 
+ 8148639 ^3 69702672 2576|125 ^ ^ 343438O8 

f^ = x^- 8x2 + 64, 7 = -2.450 - 1.414i 

(-^ _ 4027 _ 157 ^/a , 17,- /o 37 -/o 

(5 + 2V2V3, 3 + V6,9 + 3^/6) 


33/12 


y2+xy = x3 (-3-^)x2 (1^1 27^)^+ 369 63^ 

55 5 /qq 
^ ~ 576 288^'*'' 

(^7+vgT73 9+733 g + ^^ 


49/56 


y'^ + xy + y = - 874 - 171x 

^ 55 

- ~64 
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p = 2cos(f) 

(3p2 + 2p - 1, 3/92 + 2/9 - 1, /)2 + p) 


81/1 


+ xy + y = — — 95x — 697 

r< 15 

- ~64 

= x3 - 3x - 1, /9 = 1.879 

(p2+p+l,(p+l)2,(p+l)2) 


148/li 


y2 = ^3 _ (464^2 _ 32Q^ _ 1490) x'^ + x 
f^ = x^ -x'^ - ?>x + \, 7 = 2.170 

ri 1363+2927-4247^ 
- 16 

(7^+7,7^ + 27 + 1,7^+7) 


148/lii 


y2 = ^3 _ (464^2 _ 32Q^ _ 1490) x'^ + x 
f^ = x^ -x'^ -2,x+l, 7 = 0.311 

ri 1363+2927-4247^ 
- 16 

(-1272 + 87 + 40, -7^ + 7 + 4, -137^ + 97 + 42) 


148/liii 


y2 = 3.3 _ (464^2 _ 320^ _ 1490) x"^ + x 

= - x2 - 3x + 1, 7 = 1.481 

ri 1363+2927-42472 
- 16 

(7^ - 37 + 2, 72 - 27 + 1, 72 - 37 + 2) 


229/8i 


y2 = 3,3 _ (66372 - 2197 - 2485) x2 - (-3077872 + 132277 + 109691) x 

_ 205 , 105^ _ 105^2 
"-^"2^16' 47 

f^ = x^-Ax-l, 7 = 2.115 

(7 + 2, 872 + I67 + 4, 872 + 177 + 4) 


229/8ii 


y2 = ^3 _ (66372 - 2197 - 2485) x2 - (-3077872 + 132277 + 109691) x 
f^ = x^-Ax-l, 7 = -0.254 

^ _ 205 , 105 105 2 
'-^~2^16' 47 

(-372 + 7 + 13, -72 + 5, -472 + 7 + 16) 


229/8iii 


y = X — (0007 — 2197 — z4oo ) x — 1 — oU M 07 + loll /7 + iuybyi 1 x 
= - 4x - 1, 7 = -1.860 

_ 205 , 105 _ 105 2 
~ 2 ^ 16 ' 4 ' 

(72-27,72-27+1,72-27) 


725/16i 


2/2 + + 7y = + x2 - (-4477 + 4152) x - 85II67 + 59004 
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/7 
fp 



1, 7 = 1.618 



ri _ 205-3007 
- 64 



3.T + X + 1 , p 
2p^ + 5p^ - p - l,-2p3 + 5p2 



X — X 
^3 I c;^2 



-1.355 

p - 1, + 3/)2 + p) 



+ xy + jy = + x'^ - (-4477 + 4152) x - 85II67 + 59004 
2-x-l, 7 = -0.618 



^ _ 205-3007 

^ - 64 



X* - x"^ - 3x^ + X + 1, p = -0.477 



/t 

(9/9^ - 13/9^ - 21p + 19, p3 - 2p2 -2/9 + 4, 9/9^ - 13/9^ - 21p + 19) 
y2 + 1/15 (75 + 74 + 17^3 ^ 14^2 ^ ^ 

+ 462io75 (^^ + 21177676 + 4715997^ + 1829857^' + 32511857^ + 82909687^ 
+91516537 + 7962897) y - x^ - (7^ + 5203817^ + 1629947^ - 12562O7'' 

-19951007^ + 39704987^ - 31925477 + 5185452) x^ 
-TI^ (-126928637'^ + 864287877^ - 1641160677^ + 5181007157^ 
-9674266907=^ + 37575046467^ - 38928222547 + 10486471269) 



I 10639239397 -,? _ 26847806027 ^,6 , 17109768152 5 _ 4243037037.4 595858143338 ^,3 
925815 ' 925815 ' 185163 ' 28055 ' 925815 ' 



316402671731 2 , 2018798518646 
308605 ' ' 925815 ' 



2030456532221 
925815 



^ _ 2129 „7 

^ 1234420 " 29626080 



42497 „6 _ 50863 „5 
" 7406520 " 



746393 3 _ 9461581 2 _ 11911217 _ 4076357 

11850432" 59252160" 59252160" 9875360 

/^ = ,t''-x3-4x2 + 4x + 1, 7 = -1.956 
(-p3 + p2 + p + l,72-7,p2_2p+l) 



a73_ 4 
16833" 



Table 2: ramification index 3 



1/15 


+ xy + y = x^ + x^ — 
ri 55 

^ - ~54 

(5,16,20) 


135x - 660 


1/10 


+ xy + y = x^ + 26 — 


19x 


(-y 95 
^ ~ 432 

(10,6,15) 




l/6i 


y^ + xy + y = x^ + x^ — 

_ 67 
^ ~ 432 


104x + 101 
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(8,7, 14) 


l/6ii 
] 


y'^ = x^ -x"^ + IQx - 180 
C = 1/27 

(8,8,9) 


5/9 


y2 + (^1+V5j y = x^+ (if) x" (495+165V5) ^ 4125+1683V5 


^ 216 

(V5 + 3,7 + 3V5, 9 + 4v5) 


5/5 


y2 + + y = + - llOx - 880 

(~i 35 

^ ~ 108 

/ K 1 o. /i!^ c; 1 o. 7+3\/5 ^ 
l0 + ZVt),0 + ^V<3) 2 / 


8/9 


y2 = - (2^/2 + 4) x^ - (154^2 + 23l) x - 1064^2 - 1520 


n 20+10^ 
- 27 

(4\/2 + 6, 4\/2 + 6,3 + 2\/2) 


12/3 


y'^ + VSy = x^- (- V3 + 1) x^ 


c = o 

(4 + 2^3,4 + 2^3,7 + 4^3) 


13/3i 


+ xy + y = x^ + x^ - (-4957 - 637) x + 92617 + 12053 


fj = x^-x-3, 7 = 2.302 

/-< 35 

108 

(4 + V13, ii±f^,4 + V13j 


13/3ii 


y2 + xy + y = x^ + x^ - (-4957 - 637) x + 92617 + 12053 


/-y = x^ - X - 3, 7 = -1.303 

OO 

^ ~ 108 

/ ^^r. . ^ /TTT 47+13^/T3^ 

/ o+Vi.i^22 + 6V13, ii±J^vid j 


17/36 


y2 + xy + 7y = x3 - 7x2 - (196947 + 30770) x - 21455377 - 3350412 
^ _ 1+717 

' ~ 2 






^ 6545+1540aA7 
~ 432 




(5+l/^,13 + 3v^, 29+7717) 


21/3 


y2 _ ^3 ( 67''+776-3075 +474-10873+9072+427-75^ ^2 

/ 77''-4l76-2075-12674+18073-10872+5l7-63\ 

192 J ^ 
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ri 75-427-9072+3673-474+307^-77''+67^ 
^ - 144 

= + 3x6 _^ Ylx'^ -9x2 + 9, 7 = 0.770 + 0.445 
10 + 2^/7^/3, 23+5^) 


28/18 


y2 + (\/7 + l) y = x^ - (\/7 + l) x^ - (944^/7 + 2496) x + 25532^/7 + 67552 

fi 1295 1 185 /7 
^ 576 216 V ' 

(6 + 2^/7, 3 + V7, 8 + 3/7) 


49/1 


y2 + (7^ + l) y = x^ - (-7^ _ ^ _ i) ^2 - (6497^ + 9IO7 + 13l) x 
-2145I72 - 2132O7 + 6760 

/^ = x3-x2-2x + l, 7 = 1.802 

ri IO+4O7+4O72 
^ - 27 

p = 2cos(f) 

(4p2 + 3p - 1, 4p2 + 3p _ 1, p2 ^ p) 


81/1 


y2 + y = - 7 
C = 

p = -(2cos(f ))-' 


ramificat 

8/98 


ion index e = 4 
y"^ + xy + y = x^ - 55146 - 2731x 

^ 1575 
^ 256 

(3 + \/2, 20 + I2V2, 21 + 14\/2) 


8/7i 


y2 = 3,3 _ _ ^2 _ (32^ _ 48) X 

3+6^2 
^ - 16 

(8 + 4V2, 4 + ^/2, 10 + 6^/2) 


8/2i 


y2 + = 3.3 _ _ ^2 _ (^33^ + 51) X - I6OV2 - 227 

ri - 87 15 /0 
^ ~ 256 64 


8/2ii 


y2 = 42;3 _ (1116/^ + 147)x - (6966\/^ - 6859) 
C = L^5v^_^^_J^^ taken from |SijaJ 
(3 + 2^/2, 9 + 4\/2, 6 + 4^/2) 


8/7ii 


y2 = 3,3 _ (_4^ _ 14) - (-32^/2 - 48) X 
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n 3-6^ 
^ - 16 

(4 + 2^/2, 6 + 2^/2, 8 + 5^/2) 


8/2m 


y2 + = x3 - (V2 + l) - (-38^2 + 5l) x + 160^2 - 227 

^ ~ 256 ^ 64 V ^ 

= 5184a;4 + 59616a;3 + 171252x2 + 10404x + 248113, a = 0.1891 + 1.1341i 
(5 + 2^2, 6 + 4^2, 5 + 2^2) 


2624/4ii 


y"^ + xy + y = x^ - {1 + V2)x'^ - (-391^2 + 448) x + 4342^2 

ri 387 1 69, /o 
^ - ~256 + eiv^ 

fp = x'^- lOx^ + 19x2 - lOx + 1, p = 7.698 

(p,/9,p + 2^ + l) 


-6267 


2624/4i 


y2 + xy + y = - (l - a/2) - (391\/2 + 448) x - 4342^2 - 

ri 387 69 /9 
~ 256 64 

/p, =x''- 10x3 + 31x2- 30x + l, pi = 4.965 

/p2 = x^ - 41x^ + 473x2 - 1063x2 + 343x - 19, pa = 19.181 

(Pl,P2,P2) 


6267 


2304/2 


2/0 3 1 

— v3y = x'^ — 1 
C = 

p = 3 + V2VS + V2 + V3 
{p,P, p) 

Table 4: ramification index 5 




5/5i 


y2 _^ = - - (42177 + 2611) X - 1578167 - 97533 
/^ = x2-x-l, 7 = 1.6180 

1083+495\/5 
^ 200 

= 36x2 _^ ;^gg2^ _ 3299^ ^ = -56.132 
(^+1^,14 + 6^5, 16 + 7V5) 




5/180 


2/2 + xy + y = x^ - 2368 - 334x 

^ 651 
^ ~ 400 

(V5 + 3,9 + 3^5, 2^+1^) 


5/5ii 


y2 + = a-a _ ^^2 _ (4217^ + 2611) x - 1578167 - 97533 
= x2 - X - 1, 7 = -0.618 
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( 


■~i 1083+495\/5 
200 

(6 + 2^/5,4 + ^, i%7v^]) 




5/5iii 


= x^ + - 36a; - 140 

6 

^ - ~25 

(e + 2^5, 6 + 2^5, ^±1^) 


5/9 


I 
( 


/2 + xy + y = x3 + x^ + 35x - 28 

-< 3 

~ 100 

{ 9+3\/5 9+3^ 15+5^/5^ 
1^ 2 ' 2 ' 2 ^ 


725/25i 


-40897-'^ - 60017^ + 32287 + 1965 

= x"^ - - 3x2 + a; + 1, ^ ^ o.738 
^ -24+187-127^-1273 

- 25 

/p, =x^-3x3 + 4a;- 1, pi = 2.356 
fp^ = - 4x3 + 3x - 1, p2 = 3.811 

(/9l,P2,P2) 


725/25ii 


y2 + (7^ + 7) y = x^ - (-73 _ -^2 _^ 1) a;2 _ (135^3 ^ 325^2 _ i3g^ _^ 2) x 

-40897^ - 6OOI72 + 32287 + 1965 

fi=, 7 = 

ri -24+187-127^-127^ 

- 25 

/pi =x''-x3-3x2 + x + l, Pi = 2.095 
fp^ = - 8x3 _^ iQ^2 = 6.486 

(Pl,P2,P2) 


1125/5 
2 


y2 + (7 + 1) y = - (-73 + ^2 _^ ^2 _ ^2^3 _ 7^2 + 5^ + 1^ 
+673 - 1472 -27 + 12 

= x^-x3-4x2+4x + l, 7 = 1.338 
C = 

fp = x'^- 3x3 _ ^2 _^ 3^ _^ ^ ^ 2.956 
{p, P, p) 



Table 5: ramification index 6 
y2 + xy + (1 - v/3) y = x3 - (V3 + 1) x^ - (836 - 405^3) x - 4739^3 + 7704 
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n 53-387^3 
^ - 54 

(3 + ^/3, 14 + 6^/3, 15 + 8^/3) 


12/66ii 
49/91i 


7/2 + + (V3 + l) 2/ = x3 - (l - V3) x2 - (405^3 + 836) x + 4739^3 + 7704 

ri 53+3877 
^ - 54 

(6 + 2^3, 5 + V3, 9 + 4^3) 

Table 6: ramification index 7 

y2 + + ay = + a;2 - (108257^ - 244367 + 8746) x 
-9953927^ + 22354067 - 797729 

C _ 815 , 495^ 30^2 
^ - 196 + 49 ^ 1^ 

= x3 - x2 - 2a; + 1, 7 = 0.445 
p = 2cos(f) 

+ 1, 16p2 + Vlp - 8, 17p2 + I3p - 9) 


49/91ii 


y2 + + ay = + - (108257^ - 244367 + 8746) x 
-99539272 + 22354O67 - 797729 

C- _815 , 495 30 2 
^ ~ 196 ^ 49 ' 7 ' 

j^ = x^ -x^ -1x + \, 7 = -1.247 
/9 = 2cos(f) 

(p^ + + 3p _ 2, 5p2 + 3/9-2) 


49/91iii 


y2 + xy + ay = + - (108257^ - 244367 + 8746) x 
-99539272 + 22354O67 - 797729 

C_ _815 , 495 30 2 
^ 196 ^ 49 ' 7 ' 

= x^ - x2 - 2x + 1, 7 = 1.802 
/9 = 2cos(f) 


49/1 


y2 + xy + y = x^ — 70 — 36x 

fi 55 
^ 196 

p = 2cos(f) 

(2p2,2p2 + 2p, 4p2 + 3p_2) 
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Table 7: ramification index 9 



81/51i 


y"^ = x^ - (-4467^ - 8367 - 214) - (-37592172 

_ 1309 3206 1529-^2 

^'-tO z^o z^o 

= - 3x - 1, 7 = 1.879 
(p2 + 1, 4p2 + 8p + 4, 5/^2 + 9p + 3) 


- 7064017 - 


199989) X 


81/51ii 


y2 = ^3 _ (_446^2 _ 83g^ _ 214) 3^2 _ (-375921^2 

1309 3206 1529 2 
243 243 ' 243 ' 


- 7064017 - 


199989) X 










= - 3x - 1, 7 = 1.879 








/> = 2cos(f) 








(p2 + 2p + 1, p2 + 2p + 1, p2 + p + 1) 






81/51iii 


y2 = 3.3 _ (_446^2 _ g3g^ _ 214) a;2 _ (_37592i^2 


- 7064017 - 


199989) X 


^ 1309 3206 1529 2 
^ 243 243 ' 243 ^ 

= a;3 - 3x - 1, 7 = -0.607 + 1.450i 
p = 2cos(f) 

(p2 + 2p + 2, p2 + 2p + 1, p2 + 2p + 2) 







Table 8: ramification index 11 



14641/1 



+ y = x^ 



x^ - lOx - 20 

n — 14 
- "m 

p = 2cos(f^) 

(p3_2p,p3_2p,p2_l) 



The correctness of tlie suggested accessory parameters is proven whenever the Fuchsian 
group is commensurable with a triangle group, i.e. the uniformizing differential equation is a 
pullback of a hypergeometric differential equation, these 25 cases can be found in |Sija| . The 
cases 1/15 with e = 3 can be found in Krammers article |Kra96] . the case 1/6? was found 
by Elkies in [Elk98] . Whenever the quaternion algebra is defined over the rational numbers, 
Reiter determined the accessory parameter in |Rei09j . In principal it should be possible to 
adapt Krammers to prove at least a few more cases. 
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